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THE  APPLICATION  OF  GREEN’S  THEOREM  TO  THE  SOLUTION 
OF  BOUNDARY-VALUE  PROBLEMS  IN  LINEARIZED 
SUPERSONIC  WING  THEORY1 
By  Max.  A.  Heaslet  and  Harvard  Lomax 


SUMMARY 

Following  the  introduction  of  the  linearized  partial  differenr- 
tial  equation  for  nonsteady  three-dimensional  compressible  flow, 
general  methods  of  solution  are  given  for  the  two-  and  three- 
dimensional  steady— state  and  two-dimensional  unsteady— state  equations . 
It  is  also  pointed  out  that,  in  the  absence  of  thickness  effects, 
linear  theory  yields  solutions  consistent  with  the  assumptions  made 
when  applied  to  lifting— surface  problems  for  swept— back  plan  forms 
at  sonic  speeds.  The  solutions  of  the  particular  equations  are 
determined  in  all  cases  by  means  of  Green’s  theorem  and  thus  depend 
on  the  use  of  Green’s  equivalent  layer  of  sources,  sinks,  and 
doublets.  Improper  Integrals  in  the  supersonic  theory  are  treated 
by  means  of  Hadamard’s  ’’finite  part"  technique. 

Four  applications  of  the  general  solutions  are  given:  First, 
the  angle-of— attack  load  distribution  for  a  supersonic,  yawed, 
triangular  plate  with  subsonic  leading  edges  is  determined.  Second, 
downwash  is  calculated  along  the  center  line  in  the  plane  of  the 
unyawed  triangular  wing.  Third,  the  growth  of  load  distribution  is 
presented  for  subsonic  and  supersonic  two-dimensional  flat  plates 
either  starting  from  rest  at  a  uniform  velocity  or  experiencing  an 
abrupt  angle-of— attack  change.  The  transient  effects  on  lift-curve 
slope  are  then  calculated.  Finally,  the  load  distribution  and  lift- 
curve  slope  of  a  specific  swept— back  lifting  surface  are  determined 
at  a  free— stream  Mach  number  of  one. 


■•■Presented  by  Dr.  Heaslet  at  the  VII  International  Congress  of 
Applied  Mechanics,  September  19  W,  London. 
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INTRODUCTION 

If  the  effects  of  viscosity  are  assumed  small  and  shock-free 
compressible  flow  is  considered,  the  velocity  field  about  a  two—  or 
three-dimensional  body  placed  in  a  uniform  free  stream  is  irrota- 
tional  and  thus  possesses  a  velocity  potential.  In  the  determination 
of  the  pressures  exerted  on  such  a  body  or  in  the  calculation  of  the 
induced  velocity  components,  the  theoretical  aerodynamic ist  is 
concerned  essentially  with  finding  the  velocity  potential  of  the  flow 
field  and,  thus,  must  determine  the  solution  of  a  second— order  non¬ 
linear  partial  differential  equation  subject  to  certain  boundary 
conditions.  The  knofln  mathematical  difficulties  that  arise  in  the 
treatment  of  such  a  problem  make  it  expedient  to  resort  to  simplify¬ 
ing  assumptions.  In  applied  aerodynamics,  however,  efficiency  of 
flight  at  high  speeds  has  focused  attention  on  bodies  inducing 
relatively  small  velocities  throughout  the  field  of  flow  and,  as  a 
consequence,  the  demands  of  engineering  furnish  a  guide  for  the 
mathematical  simplification  of  the  theory.  The  so-called  linearized 
theory  of  compressible  flow  was  developed  to  solve  such  problems  and, 
although  considerable  work  of  a  more  precise  nature  has  been  presented 
in  two  dimensions,  a  large  amount  of  investigation  in  unsteady  flight 
f  and  in  three-dimensional  wing  theory  remains  to  be  completed  within 
the  framework  of  the  simplifying  conditions. 

The  present  paper  is  restricted  to  a  discussion  of  wing  theory 
subject  to  the  assumptions  of  linearized  compressible  flow.  It 
therefore  employs  solutions  of  Laplace* s  equation  and  the  wave 
equation  for  cases  where  the  boundary  conditions  are  specified  in 
the  plane  of  the  wing.  Attention  will  be  directed  primarily  to  the 
analysis  of  steady— state  conditions  although  an  equivalence  will  be 
established  between  the  two-dimensional  differential  equation  con¬ 
taining  the  time  variable  and  the  equation  applying  to  three- 
dimensional  supersonic  wing  theory.  Solutions  in  all  cases  will  be 
obtained  through  the  use  of  Green*  s  theorem  and  the  resultant  concept 
of  Green* s  equivalent  layer  of  sources,  sinks,  and  doublets.  The 
correspondence  between  the  theoretical  development  for  subsonic  and 
supersonic  speeds  is  particularly  useful  since  experience  related  to 
analysis  in  either  flight  regime  is  more  readily  transferred. 

In  view  of  the  widespread  use  of  sources,  sinks,  and  doublets 
in  low-speed  studies  and  the  fact  that  the  earlier  applications  to 
supersonic  wing  theory  by  Prandtl  (reference  l)  and  Schlichting 
(reference  2)  corresponded  to  the  use  of  Green* s  equivalent  layer, 
it  is  notable  that  later  emphasis  has  shifted  to  other  methods  of 
solution.  Sources  alone  were  used  by  Puckett  (reference  3)  to  create 
symmetrical  nonlifting  wings  and  were  also  applied  to  the  study  of 
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lifting  triangular  wings  with  supersonic  leading  edges ,  "but  the  use 
of  source ,  sink,  and  doublet  sheets  has  not  been  as  extensive  as 

¥  might  have. been  expected.  This  anomaly  is  even  more  apparent  in  view 

of  the  vast  mathematical  and  physical  literature  centering  around 
the  use  of  Green1  s  theorem.  One  possible  explanation  may  stem  from 
the  fact  that  the  interest  of  the  mathematician  and  physicist  in  the 
wave  equation  has  arisen  in  connection  with  problems  in  acoustics , 

*  optics,  and  vibrating  membranes.  Such  problems  introduce  boundary 
conditions  of  the  Cauchy  type,  that  is,  initial  conditions  need  to 
be  known  both  for  the  unknown  function  and  its  rate  of  change.  The 
supporting  surface  for  such  boundary  conditions  cuts  the  character¬ 
istic  cone  of  an  arbitrary  point  in  a  closed  curve  and  has  been 
called  by  Hadamard  (reference  4)  a  duly  inclined  surface.  In  aero¬ 
dynamics  the  supporting  surface  is  nonduly  inclined  and  cuts  the 
characteristic  surface  or  Mach  cone  along  the  arc  of  a  hyperbola  and, 
as  a  result,  the  problem  is  no  longer  of  the  Cauchy  type  and  the 
analysis  becomes  similar  to  that  used  in  subsonic  theory  in  the  solu¬ 
tion  of  Laplace1 s  equation.  Prior  to  the  interest  of  the  theoretical 
aerodynamic ist  in  supersonic  wing  theory,  it  appears  that  little 
attention  in  application  was  paid  to  this  type  of  solution. 

w  ^  The  material  presented  here  is  divided  into  two  main  divisions: 

Analysis  and  Applications.  In  the  first  part  of  the  Analysis 
division,  the  linearized  differential  equation  for  nonsteady  com¬ 
pressible  flow  is  given  together  with  the  underlying  assumptions  made. 
Specific  forms  of  this  equation  for  two—  and  three-dimensional  steady 
states  and  two-dimensional  unsteady  states  are  then  considered.  It 
is  pointed  out  in  particular  that  for  swept— back  lifting  surfaces 
linearized  theory  yields  consistent  solutions  at  a  free— stream  Mach 
number  of  one  although  the  analysis  of  arbitrary  thickness  distribu¬ 
tions  is  not  possible.  Following  the  various  equations.  Green1  s 
theorem  is  applied  to  find,  in  terms  of  the  known  boundary  conditions, 
the  desired  solution  by  means  of  source  and  doublet  distributions. 

Applications  of  the  general  methods  are  confined  to  four  problems . 
As  an  example  of  the  manner  in  which  angle-of— attack  load  distributions 
are  determined  for  a  lifting  flat  plate,  the  case  of  a  yawed  triangular 
wing  with  subsonic  leading  edges  is  solved.  Doublet  distributions  are 
then  applied  in  the  second  problem  to  the  calculation  of  downwash 
behind  the  same  wing  in  an  unyawed  position.  Third,  the  growth  of 
load  distribution  with  time  is  derived  for  a  supersonic  two-dimensional 
flat  plate  either  experiencing  a  sudden  sinking  motion  or  starting 
from  rest  at  a  uniform  velocity.  Such  distributions  are  of  value  in 
the  calculation  of  indicial  lift  functions  and  can  be  used,  together 
with  Duhamel1s  integral,  in  the  study  of  certain  dynamic  maneuvers. 

*  The  final  application  considers  at  a  Mach  number  equal  to  one  the  case 
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of  a  swept-back  lifting  surface  with  tips  normal  to  the  free-stream 
direction. 

All  important  symbols  are  defined  in  the  appendix. 


ANALYSIS 

The  Partial  Differential  Equations 

Basic  differential  equation.—  Consider  an  aerodynamic  body 
flying  at  an  arbitrary  Mach  number  Mq  in  air  initially  at  rest. 

If  a  Cartesian  coordinate  system  x,  y,  z  is  fixed  relative  to  the 
body,  the  body  may  then  be  assumed  stationary  and  situated  in  a  free 
stream  with  the  same  Mach  number.  If  the  free— stream  velocity  vector 
is  parallel  to  and  in  the  direction  of  the  position  x  axis  and  if  q> 
denotes  the  perturbation  velocity  potential  for  isentropic  flow,  the 
linearized  partial  differential  equation  for  cp  may  be  written  in 
the  form 


(Mo2—1)  ^XX  “  Tyy  —  Tzz  +  Ttt  +  — ^  <Pxt  =  0  (l) 

ao2  ao 

where  a0  is  the  velocity  of  sound  in  the  free  stream  and  t  denotes 
time. 


The  assumptions  underlying  the  derivation  of  equation  (l)  have 
been  stated  in  numerous  places  but  are  not  always  obviously  compatible. 

It  is  assumed  here  that  the  ratios  are  small  compared  to 

,  *0  V0  v0 

one,  where  u,  v,  w  are  induced  velocity  components  and  V0  is  the 
velocity  of  the  free  stream;  moreover. 


7~1  w  2  f 211  U2+V2+W2 


«1 


and,  finally,  the  velocity  gradients  at  a  given  point  of  the  flow 
field  are  all  of  similar  magnitude. 


Special  cases.—  The  particular  forms  of  equation  (l)  to  he 
considered  are  given  i«n  table  I.  In  the  steady-state  equations  the 
original  independent  variables  are  retained;  the  two-dimensional 
unsteady-state  equation  has,  however,  new  variables  defined  by  the 
relations 


NACA  TN  No.  1J6J 


5 


y 


x*  =  x  -  aoMot,  z*  =  z,  t*  =  a0t  (2) 


and  consequently  the  "boundary  conditions  for  any  particular  example 
will  "be  subject  to  the  same  transformation.  In  all  the  equations 
the  constraints  imposed  "by  the  linearization  permit,  for  problems  in 
wing  theory,  the  boundary  conditions  to  be  specified  in  the  plane  of 
the  wing.  This  plane  shall  arbitrarily  be  taken  to  be  z=0. 


TABLE  I.-  LINEARIZED  PARTIAL  DIFFERENTIAL  EQUATIONS  OF 

COMPRESSIBLE  FLOW 


Steady  State 

Two  J 

(l“M02)  9xx  +  ^zz  "  Oj 

Mo<l, 

(A) 

Dimensions 

(Mo2—1)  9xx  —  Tzz  =  0, 

Mo>  1, 

(B) 

"(i-Mo2)  cp^  +  qiyy  +  cpzz 

=  0,  M0<1, 

(c) 

Three  h 

(Mo2— l)  Txx  “  cPyy  ~  ^zz 

=  0,  Mq>1, 

(B) 

Dimensions 

<-cpyy  +  ^zz  =  0 

Mq  =  1, 

(E) 

Unsteady  State 

Two 

Dimensions 

Tt*t*  ~  Txfxt  —  T  z f  z  *  - 

0 

(F) 

The  Mach  number  range  for  which  the  equations  are  valid  cannot 
be  prescribed  a  priori  since  induced  velocities  are  functions  of  wing 
geometry  and  angle  of  attack.  We  can  say,  however,  that  for  certain 
configurations  at  small  angles  of  attack  the  equations  and  their 
solutions  are  consistent  with  the  assumptions.  In  particular,  three- 
dimensional  lifting  surfaces  with  sufficient  sweepback  yield  solutions 
of  this  class  at  Mo=l.  The  differential  equation  shows  that  in  this 
case  the  boundary  conditions  need  only  be  specified  along  strips  in 
the  transverse  direction.  The  surfaces  of  the  Mach  cones  also  are 
normal  to  the  free-stream  direction  so  that  any  disturbance  point 
makes  itself  felt  at  all  points  not  upstream  of  it.  Since  for  these 
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lifting  surfaces  the  disturbances  do  not  become  excessive  at  Mo=l, 

we  have  a  specific  kind  of  lateral  strip  theory  that  yields  formal 

solutions  compatible  with  the  assumptions  made.  * 


Boundary  Conditions 

Steady  state.—  The  boundary  conditions  are  given  in  the  z=0 
plane  and  in  the  case  of  two-dimensional  theory  the  wing  is  assumed 
to  extend  infinitely,  parallel  to  the  y  axis.  As  a  convenience  in 
notation,  two  subscripts  will  be  introduced:  the  first,  u,  denotes 
conditions  on  the  upper  surface  of  the  wing,  that  is,  the  limit  of 
the  function  as  z  approaches  zero  through  positive  values;  the 
second,  Z,  denotes  conditions  on  the  lower  surface  of  the  wing,  that 
is,  the  limit  of  the  function  as  z  approaches  zero  through  nega¬ 
tive  values. 

Four  types  of  boundary  conditions  arise  in  practice: 

1.  Symmetrical  nonlifting  wing  (boundary— value  problem  of  the 

first  kind) The  conditions  wu  =  vj  =  0  hold  over  all 
of  the  xy  plane  except  for  the  region  occupied  by  the  wing. 
On  the  wing,  the  relations  2wu  =  —  2w^  =  Aw0  =  f(x,y) 
are  given,  the  function  f(x,y)  being  determined  by  the 
geometry  of  the  configuration.  Over  all  of  the  xy  plane, 
Au0  =  uu  —  u^  =0  applies. 

2.  Lifting  surface  with  specified  loading  (boundary  value 

problem  of  the  first  kind) .—  The  condition 
Au0  =  uu  —  vli  =  0  holds  over  the  xy  plane  except  for  the 
region  occupied  by  the  wing.  On  the  wing,  the  relations 
2uu  =  -  2u i  =  Au0  =  f(x,y)  are  given,  the  function  f(x,y) 
being  determined  by  the  specified  loading.  Over  all  of 
the  xy  plane,  Aw0  =  0  applies. 

3.  Lifting  surface  with  specified  camber  and  angle  of  attack 

(boundary— value  problem  of  the  second  kind)  The 
condition  Au0  =  0  holds  over  the  xy  plane  except  for 
the  region  occupied  by  the  wing.  On  the  wing  the 
relation  w  =  f(x,y)  is  given,  the  function  f(x,y) 
being  determined  by  the  given  camber,  twist,  and  angle 
of  incidence.  Over  all  of  the  xy  plane,  Aw0  =  0  applies. 

4.  Symmetrical  wing  with  specified  pressure  distribution 

(boundary— value  problem  of  the  second  kind) .—  The 
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condition  Awq  =  0  holds  over  all  of  the  xy  plane  except 
for  the  region  occupied  hy  the  wing.  On  the  wing  the 
relation  u  =  f(x,y)  is  given,  the  function  f(x,y)  being 
determined  hy  the  specified  pressure  distribution.  Over 
all  of  the  xy  plane,  Au0  =  0  applies. 


In  all  cases,  induced  velocity  u  is  related  to  pressure  coef¬ 
ficient  Cp  by  the  relation 


Unsteady  state.—  The  steady— state  boundary  conditions  have  been 
given  in  the  most  general  terms  possible.  The 'unsteady-state  condi¬ 
tions  will  be  limited  to  a  more  restricted  type  of  problem,  namely, 
cases  wherein  the  airfoil  is  assumed  to  experience  at  t=0  either  an 
abrupt  change  in  angle  of  attack  without  pitching  or  starts  to  travel 
at  the  instantaneous  velocity  V0  and  angle  of  attack  a.  In  this 
way  the  transient  variation  of  load  distribution  and  airfoil  charac¬ 
teristics  can  be  calculated  for  unit  angle— of-attack  change.  Similar 
methods  can  treat  unit  rate  of  pitching,  or  deflection  of  aileron,  as 


well  as  the  effects  produced 
when  the  airfoil  enters  a  gust 
of  given  structure.  The  use 
of  solutions  of  such  problems 
in  connection  with  operational 
methods  is  well  known  in  applied 
mathematics.  Applications  of 
these  operational  methods  to 
aerodynamics  have  been  given  by 
R.  T.  Jones  (references  5  and  6) 
for  incompressible  fluid  theory 
and  by  Heaslet  and  Lomax  (refer¬ 
ence  7)  for  supersonic  flow. 

If  the  rectangular  coor¬ 
dinate  system  xl,zs,t11  assoc¬ 
iated  with  equation  (2)  is  con¬ 
sidered  to  be  fixed,  the  airfoil 
moves  in  the  negative  x*  direc¬ 
tion  and  the  free— stream  velocity 
is  zero.  A  simple  distortion  of 
the  time  axis  is  also  introduced 
to  simplify  the  differential 
equation.  Figures  l(a)  and  1(b) 
aid  in  the  visualization  of  the 


(a)  Supersonic  wing. 


(b)  Subsonic  wing. 


Figure  1.-  Diagram  for  use  in  deter¬ 
mining  boundary  conditions  in  two- 
dimensional  unsteady  motion. 
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problems  involved.  The  airfoil  section  is  assumed  to  lie  initially 

on  the  x’  axis  with  leading  edge  at  the  origin  and  trailing  edge  at 

the  point  xJ  =  c0.  As  time  progresses  the  airfoil  sweeps  across  a  ^ 

portion  of  the  x’t’  plane,  the  leading  edge  traversing  the  line 

x’  =  —  M0t’  while  the  trace  of  the  trailing  edge  is  the  line 

x*  =  Co  —  Mot* .  The  region  bounded  by  these  lines  and  the  line 

t’  =  0  is  that  swept  by  the  airfoil.  The  characteristic  cones  of 

the  differential  equation  cut  the  x’t’  plane  along  lines  inclined 

at  ±  ^5°  to  either  axis .  If  the  airfoil  experiences  an  angle— of— 

attack  change  a  without  pitching,  the  "area"  swept  over  by  the  wing 

must  yield  w  =  —  V0a.  On  the  other  hand,  if  the  airfoil  enters  a 

gust  of  constant  vertical  velocity  w0,  the  region  over  which  the 

modification  of  w  is  effective  is  restricted  to  the  region  occupied 

simultaneously  by  the  airfoil  and  the  gust.  If,  for  example,  the 

edge  of  the  gust  is  fixed  along  the  t*  axis,  this  axis  will  form  the 

right-hand  boundary  of  the  region  over  which  the  change  in  boundary 

conditions  occurs.  A  statement  of  these  boundary  conditions  may  be 

put  in  the  following  form:  * 

1.  Lifting  surface  undergoing  abrupt  change  or  starting  from 

rest  with  given  velocity.—  The  condition  Au0  =  0  holds 
over  all  of  the  x’t’  plane  except  for  the  region  swept 
across  by  the  airfoil.  In  this  latter  region,  the  rela¬ 
tion  w  =  f(x* ,t’)  is  given,  the  function  f(x’,t’)  being 
determined  by  the  modification  in  airfoil  angle  of  attack, 
pitching  velocity,  aileron  deflection,  or  by  the  gust 
structure.  Over  all  of  the  x’t’  plane,  Aw0  =  0  applies. 

The  expression  for  pressure  coefficient  is 


c  .J*.  V  ^  ) 

Vc/-'  V  Ot  0  ox  J 

=  -2  dcp 
VqM0  8t*’ 


Solution  of  Boundary— Value  Problems 

General  treatment.-  The  use  of  Green’s  theorem  in  the  solution 
of  second— order  partial  differential  equations  leads  one  to  the  con¬ 
sideration  of  certain  particular  solutions  of  the  given  equations.  * 
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Because  of  the  physical  importance  as  well  as  the  mathematical  appli¬ 
cability,  attention  has  been  centered  on  the  use  of  a  so-called 
fundamental  solution  or  source  potential.  Thus,  in  the  subsonic 
case,  the  potential  at  the  point  x, z  of  a  unit  source  located  at 
the  point  |,  £  and  applicable  to  equation  (A),  table  I,  is  the 
logarithmic  function 


<p(x,z)  =  7  In  rA  =  ^  In  <v/(x— |)2+p2( z-£)2  (3) 

2jt  a  2jt 


* 


while  for  equation  (C)  the  potential  at  x,y,z  of  the  unit  source 
at  is 


cp(x,y,z) 


-1  =  _ =1 _ 

kitrc  4jt  a/C x-i  )  2+P2  [  ( y-rj )  +  ( z-0  2  ] 


(4) 


Here  and  elsewhere  we  have  p2  =  1 1-Mq2  |  where  the  bars  indicate 
that  absolute  values  are  to  be  taken. 


The  application  of  these  potential  functions  to  the  solution  of 
boundary— value  problems  in  subsonic  linearized  flow  is  well  known. 
Supersonic  theory,  however,  introduces  added  complications  when 
the  fundamental  solutions  are  considered  and,  although  methods  have 
been  established,  the  mathematical  techniques  are  of  comparatively 
recent  origin.  The  principal  difficulty  lies  in  the  integration  of 
higher-ordered  singularities  that  appear  in  the  three-dimensional 
analysis.  Hadamard  (reference  4)  resolved  these  difficulties  and 
thus  avoided  the  more  specialized  approach  of  Volterra  (reference  8). 
It  would  appear,  however,  that  a  more  direct  method  of  derivation 
stems  from  Marcel  Riesz*s  use  of  fractional  integration.  (See,  in 
this  connection,  references  9  and  10.)  The  oddness  or  evenness  of 
the  number  of  dimensions  still  involves  considerable  differences 
but  the  final  solutions  are  easily  applied. 

In  two-dimensional  supersonic  flow,  the  potential  at  the  point 
x, z  of  a  unit  source  located  at  the  point  |,  £  is  defined  as 
follows: 


cp(x,z)  =  0  for  (x-|)2<p2(z-£)2 
<p(x,z)  =  -  T~  for  (x-?)2^  p2(z-^)2 


W 


(5) 
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In  three  dimensions,  the  source  potential  is 


<P(x,y,z) 


-1  =  -1 
2orrD  2jt  a/  (x— |)2— (32[  (y— ti)2+(  z— ^  )2  ] 


(6) 


ir 


at  all  points  for  which  the  radical  is  real  and  is  zero  elsewhere. 

These  functions  are  directly  applicable  to  equations  (B)  and  (D) 
of  table  I.  Equations  (E)  and  (F)  are,  of  course,  special  mathe¬ 
matical  cases  of  equations  (A)  and  (D)  for  which  Mo  is  0  andV^? 
respectively. 

By  means  of  the  various  source  potentials,  it  is  now  possible 
to  present  solutions  of  the  differential  equations  in  terms  of  the 
prescribed  boundary  conditions.  These  conditions  are  assumed  to  be 
given  in  the  z=0  plane  and  subscripts  u  and  l  shall  again 
denote  conditions  at  z=0+  and  z=0-,  respectively.  The  general 
solutions  appear  as  follows: 

Equation  (A),  Mq<  1, 


<p(x,z) 


Equation  (B),  Mq>1 


cp(x,z) 


i  rx_pz  ScPu 

Pj  dz 

—  00 


for  z  >  0 


<S| 


for  z  <  0 


(8) 


Equation  (C),  Mq<1 


<P(x,y,z) 


Sep 

Sz 


0^=0 


d?dT] 

(9) 


<T 


w 
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Equation  (D),  Mq>1 


<P(x,y,z) 


'l_\  /ZVa 

^4=0  \  Sz  Sz ) 


(10) 


In  the  last  equation,  the  range  of  integration  is  confined  to  that 
portion  of  the  z=0  plane  that  lies  within  the  Mach  forecone  of  the 
point  x,y,  z,  that  is,  within  the  half— portion  of  the  right  circular 
cone 


(x-|)2-02[(y-^)2+(z-t)2]  =  0 


lying  upstream  of  the  point  x,y,  z.  The  semi  vertex  angle  p  of  this 
cone  is  the  Mach  angle  and  is  given  "by  the  relation 


p  =  arc  sin  —  =  arc  cot  6 
M0 

The  symbol  [  was  introduced  by  Hadamard  and  denotes  the 

"finite  part"  of  the  integral.  As  in  the  case  of  Cauchy1  s  principal 
value,  an  improper  integral  is  reduced  by  a  prescribed  technique  to 
a  finite  and  unique  value .  By  definition  ( see  also  reference  11) , 


nX°  A(x)dx  =  n  X°  A(x)-A(xq)  _  2A(xq) 
a  (x0-x)3/2  Ja  (x0-x)3/2  (xQ-a)1/2 


(ID 


In  the  two-dimensional  supersonic  case,  the  solution  for  the 
velocity  potential  is  expressed  as  the  integral  of  a  distribution  of 
source  potentials.  In  all  other  cases,  both  the  source  potential 
and  its  derivative  appear  in  the  integrand,  this  latter  expression 
being  identified  with  the  doublet  potential. 

Nonlifting  case  (symmetrical  wing,  boundary— value  problem  of  the 
first  kind.—  Equations  (7)  *  (8),  (9) /and  (10)  are  applicable  directly 
to  the  calculation  of  the  potential  function  corresponding  to  a 
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symmetrical  wing.  The  relation  cpu  =  Tz  follows  from  the  condition 
uu  •=  uz .  Moreover,  if  — denotes  the  local  slope  of  the  upper 
surface  of  the  wing, 


2_  <*Pji  =  _  i_  =  /dZjA 
VG  5  z  V0  Sz  \  dx / 


* 


1 


and  the  solutions  of  the  various  equations  are  expressed  in  terms  of 
source  distributions  alone. 

For  example,  equation  (10)  becomes 


<p(x,y,z) 


dzu  d£dr] 

^  </  (x— |)2— 02[(y-T))2+z2] 


(12) 


where  the  finite  part  sign  is  dropped  since  the  integral  is  proper. 
This  equation  was  given  by  Puckett  in  reference  3* 

The  pressure  coefficient  on  the  surface  of  the  wing  is 


cr  '  "  %  dCLo 

Lifting  case  (boundary— value  problem  of  the  first  kind)  From 
the  condition  wu  =  w^  we  have  Aw0  =  0  and  the  integrands  in 
equations  (7)*  (9)>  and  (10)  are  expressed  solely  in  terms  of 
doublet  distributions,  while  equation  (8)  yields  the  result  that 
conditions  on  either  side  of  the  wing  have  no  effect  on  the  other 
side. 


Taking  equation  (10)  as  an  example,  the  solution  under  the  pre¬ 
scribed  conditions  is 


cp(x,y,z)  =  -  ~ 


f  J 


|3gz  dgdri 


(x-5)2-P2[(y-^)2+z2]  a/2 


(13) 
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A  more  direct  evaluation  of  perturbation  velocity  u  can  be  obtained 
from  the  alternate  expression 


(i*0 


Similar  expressions  exist  for  equations  (7)  an(l  (9) • 

Lifting  case  (boundary— value  problem  of  the  second  kind) This 
type  of  boundary  condition  cannot  be  solved  directly  by  means  of  the 
formulas  which  have  been  presented  but  resolves  always  into  the 
required  solution  of  an  integral  equation.  In  three-dimensional  sub¬ 
sonic  wing  theory,  the  method  of  solution  depends  usually  on  seme 
modification  of  Prandtl’s  lifting-line  theory  although,  more  recently, 
lifting-surface  theories  by  Falkner  (reference  12)  and  Cohen  (refer¬ 
ence  13)  have  been  applied  successfully. 

In  the  case  of  three-dimensional  supersonic— wing  theory,  sources, 
sinks,  and  doublets  have  been  utilized  in  two  ways  in  the  solution 
of  lifting-surface  problems.  The  first  of  these  methods  was  given  by 
Eward  (references  ik  and  15)  and  is  particularly  powerful  when  one  of 
the  leading  edges  of  the  wing  is  of  the  supersonic  type,  that  is, 
when  the  velocity  component  of  the  free  stream  normal  to  the  edge  is 
greater  than  the  speed  of  sound.  A  second  method  of  solution  was 
presented  in  reference  1 6  for  the  important  case  of  wings  with  sub¬ 
sonic  leading  edges,  provided  the  flow  field  about  the  wing  is  of  the 
conical  type  introduced  by  Busemann  (reference  17).  The  essential 
feature  of  this  method  is  the  use  of  differential  lifting  elements 
carrying  a  constant  load  and  designed  for  use  in  conical  flow  fields. 
The  solution  consists  of  determining  the  distribution  of  loading 
over  these  elements  so  that  the  resultant  induced  vertical  velocity 
at  any  point  on  the  lifting  surface  satisfies  the  local  boundary 
condition.  When  approached  from  this  standpoint,  the  problem  again 
requires  the  solution  of  an  integral  equation  but  the  equation  is  of 
the  form 

w(x)  = 


f (xi)dxi 
Xi— x 
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and  is  thus  well  known  from  low-speed  airfoil— section  theory. 
Inversions  of  this  equation  have  been  provided  by  Allen  (refer¬ 
ence  18)  and  von  Mises  and 
Friedrichs  (reference  19) . 


Ap 

q 


M0 


4  •'  #  4 

S+a8 

y 


Figure  2  shows  the 
elemental  lifting  surface  to  be 
used.  The  sides  of  the  element 
extend  back  from  the  tip  of  the 
Mach  cone,  making  angles  5  and 
5+A6  with  the  positive  x  axis 
or  free— stream  direction.  The 
vertical  velocity  induced  at  the 
point  x,y,  0  by  the  element  will 
be  a  function  of  5,  AS,  and  ^ 

or,  changing  the  notation,  0,  A 9 


i 


Figure  2.—  Lifting  surface  element 
carrying  constant  load. 

and  od  vhere 


Q  =  p  tan  b 
6+AQ  =  3  tan  (5  +  A&) 


Denoting  the  gradient  of  this  vertical  velocity  with  respect  to  0 
8wz_o 

by  — — ,  it  follows  that 

80 


^2=0  _  lim  w(0+A0,co)-w(0,a>) 

A0-^O  A0 

where  w(0,o>)  and  w(04A0,<n)  are  the  velocities  induced  by  right- 
triangular  lifting  surfaces  with  constant  loading  and  with  vertex 
angles  equal  to  5  and  6+A&,  respectively.  The  velocities 
induced  by  the  constantly  loaded  surface  are  determined  directly 
from  equation  ( 13) -  The  results  of  these  calculations  yield  the 
following  expressions: 


For  cjd  <  0 
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^z=o  =  _  _cp_  rm  Vw 

50  2  rtV0  J_x  CDi(a>i-0)2 

and  for  o>  >  0 


(16a) 


2s=°  -  -  -£ g_  r"  aioj  (16b) 

50  2jtV0  ■A  mi(<nx-0)2 

■where 

C=M4E 

2  q 

The  term  ^  is  the  load  coefficient  and  is  equal  to  the  difference 

between  pressure  coefficients  on  the  lower  and  upper  surface  of  the 
wing, 

Ap  _  pz-pu 

T - q - %  "  CPu 


Nonlifting  case  (boundary-value  problem  of  the  second  kind).—  In 
the  previously  discussed  lifting  case,  the  induced  vertical— velocity 
field  of  a  constantly  loaded  element  was  calculated.  An  analogous 
type  of  element  can  also  be  developed  for  use  in  the  determination  of 
nonlifting  wings  with  prescribed  pressure  distributions.  It  is 
apparent  that  differential  expressions  similar  to  equations  (l6a)  and 
(l6b)  must  be  derived  which  establish  the  induced  field  of  the 
x  component  of  perturbation  velocity  for  a  conical— flow  element  with 
constant  vertical  velocity.  From  equation  (12),  the  following 
expressions  result: 

For  cd  <  0 


5Cp  _  2\  P®  .  (Didcoj. 

50  _ x  ((Ox— 0)2  V 1— CDx2 


(17a) 


I 


1 6 

and  for  m  >  6 
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3Cp  _  2X  ft  (Didmi  (17b) 

be  3tt  ^  (a)i-0)2  V1-cui2 

where  pressure  coefficient  Cp  and  surface  slope  X  are 


The  application  of  equations  (17a)  and  (17b)  to  the  determination 
of  a  thickness  distribution  supporting  a  given  pressure  distribution 
consists  of  determining  X  as  a  function  of  0  such  that  the  desired 
pressure  results.  The  essential  simplification  of  the  method  is  brought 
about  by  the  use  of  elements  which  lead  to  single  integral  equations  of 
standard  form. 


APPLICATIONS 

Yawed  Triangular  Lifting  Surface 

Consider  a  yawed,  triangular  flat  plate  with  subsonic  leading 
edges  such  as  is  indicated  in  figure  3*  Relative  to  the  x  axis  or 


free— stream  direction,  the  sides  of 


Figure  3«-  Yawed  triangular  flat- 
plate  lifting  surface  with  sub¬ 
sonic  leading  edges. 


the  triangle  make  the  angles 
60  and  b±  so  that  the  total 
vertex  angle  is  b0+bx=2A. 

The  quantities  60,  0 i,  and  u> 
are  also  introduced  where 


0O  =  P  tan  50 


0!  =  (3  tan 


4 


(JO 
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The  vertical  induced  velocity  at  any  point  on  the  wing  can  he 
found  from  equations  (l6a)  and  ( l6b ) ,  provided  the  distribution  of 
the  loading  factor  C  is  known.  Setting  C=C(0) ,  the  downwash  is 
given  by  the  expression: 


18 
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Substitution  of  equation  (22)  into  equation  (l8)  yields  the  two 
relations 


and 


where 


wz=0  -  2^  [AHi(©o>®i)  +  BH2(0q>0i)  ] 


wz=0  =  ay"  +  BH2(0x,0o)] 


Si(0i^0q)  = 


dn)i 


03i  J (l-O)!2)  (cDi—Oi)  (a>i+0o) 


and 


H2(0ij3o)  - 


0i 


dn)i 


0)l2  \/  (l~<X>i2)  (<Di+0O) 


(23) 

(24) 

(25) 


(26) 


Equations  (25)  and  (26)  may  be  integrated  by  the  standard  methods 
for  elliptic  integrals  and,  after  substituting  into  equations  (23) 
and  (24)  and  solving  for  A  and  B,  we  have 


A 


lo^O  (  /3I 

(BE*  V  0Q+01 


(27) 


and 

B  =  —  --Q— 5--°  20o0i  /  — 

|3E*  v  ©o*^"0i 

where 

G  =  1+aoQi  ~  y(l-0o2)(l-0i?) 

0Q+01 


(28) 


(29) 


and  E*  is  the  complete  elliptic  integral  of  the  second  kind  with 
modulus  v 1 ~G2 . 

The  load  distribution  over  the  wing  can  now  be  calculated  from 
equations  (22),  (27),  and  (28).  It  follows  that 
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4 


S 


ap  =  gc(0)  =  j±  /ZHZ 

q  Vo2  3E»a/  00+01 


(  0q”^i)  0+20Q0; 


(30) 


Typical  load  distributions  over  a  yawed  wing  are  shown  in  figure  k  for 
3  tan  61  =  0.6  and  for  3  tan  50  equal  to  0,  0.3>  0.6,  and  0.9. 


-1.0  -.8  -.6  -4  -.2  0  .2  .4  .6  .8  1.0 


j9  TAN  8 

Figure  4.—  Angle— of— attack  load  distribution  over  yawed  triangular 
plan  form,  3  tan  61  =  0.6.  - 
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Integration  of  equation  (30)  over  the  surface  of  the  triangular 
wing  determines  for  lift  coefficient  the  expression 

CL  ,  SM  cosAyOip  (3!) 

where  A  is  sideslip  angle  and  2A  is  the  angle  between  the  leading 
edges.  Equation  (31)  was  derived  for  wings  with  subsonic  leading 
edges  and  supersonic  trailing  edge  and  consequently  is  valid  only  for 
cases  for  which 


p  +  A<  90° 

A  +  A<  p 

A  -  A>  0  (32) 


Downwash  Behind  Triangular  Wing  in  Supersonic  Flow 


The  second  application  will  show  how  doublet  distributions  may 
he  employed  in  the  calculation  of  downwash  in  the  wake  of  an  unyawed 
triangular  wing  with  subsonic  leading  edges.  The  expression  for  the 

velocity  potential  will  be 
given  in  general  form,  but 
in  order  to  avoid  detailed 
analysis  the  value  of  down- 
wash  is  determined  only 
along  the  center  line  in 
the  plane  of  the  wing. 

A  plan  view  of  the 
wing  and  wake  is  shown  in 
figure  5.  The  load  distri¬ 
bution  over  the  wing  is 
found  from  equation  (30), 
after  setting  0o~01}  to 
be 

40o2ax  _ 


AP  =  _ 
q  E‘op  Ve02x2-p2y2 


(33) 


Figure  5-—  Triangular  wing  and  wake 
together  with  regions  used  in  cal¬ 
culating  downwash. 
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where  E’Q  has  the  modulus  a/l— 0O2 .  Moreover 


* -r  u 

J-00 


dx 


so  that  the  jump  in  potential  cj^-cp^  on  the  surface  of  the  wing  is 

<Pu-9Z  =  a/ So2*2-?2;/2  (34) 

E’oP 

and  for  points  in  the  wake  is 


<*rn  - 

&  0? 


(35) 


where  cQ  is  the  root  chord  of  the  wing.  Since  the  wing  and  its 
wake  form  a  discontinuity  surface  for  the  perturbation  velocity 
potential  and  since  for  all  points  on  this  surface 

dq?u  _  &p  i 

Bz  5z 


it  follows  from  equation  (10)  that  the  velocity  potential  at  an 
arbitrary  point  x,y,  z  is  given  by  the  relation 


9(x,y,z)  = 


lia” 


(36) 


where  t  is  that  portion  of  the  wing  and  wake  forward  of  the  Mach 
forecone  from  the  point  x,y,  z  and  (fpu- ^z)  is  given  by  equa¬ 
tions  (34)  and  (35). 

The  value  of  the  downwash  aft  of  the  wing  and  along  the  x  axis 
will  be  calculated  from  equation  (36),  thus 
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w  = 


a 

3^ 


<p(x,y,z) 


z=0 

y^=o 


« 


is  to  "be  determined.  In  carrying  out  these  calculations,  it  is 
necessary  to  consider  two  segments  of  the  x  axis  behind  the  wing; 

Region  A  extends  from  the  trailing  edge  to  the  point  where  the 
trailing  Mach  cones  from  the  tips  of  the  wing  intersect  the  x  axis 
and  thus  includes  values  of  x  satisfying  the  inequality 


c0^x^  c0(l+9o) 


Region  B  includes  values  of  x  for  which 


co(l+0o)  ^  x 


The  final  expressions  for  downwash  in  the  two  regions  are 
found,  after  some  manipulation,  to  he 


Region  A; 


JL  =  2  E2-(1-~k22)K2  +  E’o-ep 

wQ  jtE’oka  E*o 

-  2  rk2  K-E  d 

JtE«o  Jo  k2(l+0ok) 


(37) 


Region  B; 


_w  _  2Ei  +  2  K-E 

wo  rtE’o  rtE'o  J0  k+0o 


(38) 


where 


w0  induced  vertical  velocity  on  the  wing 

K  complete  elliptic  integral  of  the  first  kind 

Ki,  K2  complete  elliptic  integral  of  the  first  kind  with 
moduli  kj.,  k2,  respectively 


E 


complete  elliptic  integral  of  the  second  kind 
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Ex,  E2  complete  elliptic  integral  of  the  second  kind  with  moduli 
kx,  k2,  respectively 

k  £o£o 

X-C0 

k2  ££a 

Cq@0 

Figure  6  shows  the  variation  of  along  the  x  axis  for  various 

wo 

values  of  the  parameter  0o=3  "tan  §o-  The  asymptotic  values  at  x=oo 
are  also  indicated  and  can  “be  shown  to  agree  with  the  values  of 
downwash  at  infinity  for  a  wing  with  the  same  span  load  distribution 
in  incompressible  flow.  The  discontinuity  in  downwash  at  the  trailing 
edge  is  a  characteristic  property  of  supersonic— type  trailing  edges. 


■jr  ,  DISTANCE  IN  CHORDS 


Figure  6.-  Variation  of  downwash  on  x  axis  behind  a  triangular 
wing  plotted  as  a  function  of  distance  in  chord  lengths. 


2k 
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Denoting  downwash  at  this  point  "by  wt,  Lagerstrom  and  Graham 
(reference  20 )  have  shown  that 


vt  _  Ey Q—0Q 

wQ  E*o 


(39) 


A  more  detailed  development  of  the  results  shown  in  figure  6  has 
heen  given  in  reference  21. 


Two-Dimensional  Unsteady  Lift  Problems 

It  has  already  been  pointed  out  that,  in  the  case  of  unsteady 
motion  in  a  two-dimensional  compressible— flow  field,  the  linearized 
partial— differential  equation  for  the  perturbation  velocity  potential 
can  be  transformed  into  the  same  form  that  has  been  considered  in 
solving  steady— state  problems  in  supersonic  wing  theory.  This 
immediately  indicates  the  possibility  that  for  certain  types  of 
boundary— value  problems  in  the  unsteady  case  an  analogy  can  be  estab¬ 
lished  with  three-dimensional  lifting— surface  problems. 


As  an  example,  let  us  consider  an  airfoil  that  has  been  flying 
at  supersonic  speed  and  then  experiences  at  tf=0  an  abrupt  angle— 
of— attack  change  without  pitching.  Since  the  angle— of— attack  change 
is  as  Slimed  to  take  place  at  t*=0,  it  can  be  assumed  that  previous 
to  this  time  the  induced  velocities  of  the  wing  are  zero  and  only 
subsequent  perturbations  are  to  be  calculated.  Throughout  the  swept 
area  in  the  xfty  plane  (fig.  l)  the  vertical  induced  velocity  w^is 

constant  and  equal  to  —Vo ct.  Elsewhere  in  the  z,=0  plane  there  is 

ocp 

no  discontinuity  in  the  value  of  pressure,  that  is, 
uous  at  z 7 =0 . 


at* 


is  contin- 


Suppose  now  that  the  area  is  a  wing  plan  form  and  that  the 
free  stream  is  directed  along  the  t*  axis.  The  characteristic 
cones  of  the  unsteady  problem  become  the  Mach  cones  of  the  steady- 
state  problem,  and  the  Mach  number  of  the  free  stream  is  ^f2  since 
the  characteristic  lines  in  the  figure  are  inclined  V5°  to  the  axes. 

a© 

Moreover,  the  induced  vertical  velocity  is  t— -  and  the  perturba— 

oz'  Sep 

tion  velocity  in  the  free— stream  direction  is  ^7*  A  correspond¬ 
ence  can  thus  be  established  between  the  unsteady  problem  and  a 
three— dimensional  lifting— surface  problem. 
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As  outlined,  the  boundary-value  problem  is  one  of  the  second 
kind,  that  is,  w  is  specified  on  the  wing  and  Au=Aw=0  off  the 
wing.  In  our  particular  example,  however,  the  edges  of  the  wing  are 
of  the  supersonic  type  and  no  interaction  exists  between  the  two 
surfaces  of  the  lifting  plate  so  that  pressures  on  either  side  can  be 
calculated  by  the  methods  used  for  symmetrical  nonlifting  wings. 

Thus,  from  equation  (12),  for  z’  >0 


cp(t»,x*,z») 


_ dtt1dxt1 _ 

y  (t*— t1!)2— (x*— x’l)2— Zf2 


(to). 


and  for  all  z’ 


<p(t*,x*  ,zr)  =  —  q>(t*,x*,— : z») 


The  expressions  for  the  indicial  load  coefficient 
follows : 


4p 

q 


are  as 


Region  A  (between  lines  x'=t»,  t»=0,  and  x»=c0-MotJ) 


Ap  4a 
q  ~  Mo 


(4la) 


Region  B  (between  lines  x'  =  -t»,  x»=t,  and  x»=c0-M0t*) 


Ap  _  4a 
q  "/Mq2— 1 


1  Mnx»+t* 

—  arc  cos  — -  + 


jtMo 


Region  C  (between  lines  x*  =  — Mot’ 


and  x’  =  — t*) 
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The  growth  of  —■  with  time,  as  obtained  from  equations  (4la), 

(4lb),  and  (4lc),  is  shown  in  the  portion  of  figure  7  designed 

"supersonic.”  At  tf=0  the  loading  jumps  to  the  value  and  is 

Mq 

constant  along  the  entire  chord.  This  value  persists  throughout  the 
previously  denoted  Region  A  and  thus,  with  advancing  time,  moves 
rearward  along  the  chord,  leaving  the  trailing  edge  at  tT  = 

Over  the  forward  portion  of  the  chord  the  familiar  Ackeret  type  of 
steady— state  loading  becomes  effective,  spreading  back  from  the 

leading  edge  and  occupying  the  entire  chord  length  after  tf  =  —rr. 

c0  Mo~ 1 

Previous  to  t1  =  — — -  a  transition  region  between  the  two  types  of 

Motl 

constant  loading  exists,  and  subsequent  to  this  time  this  transition 


region  moves  aft  and  leaves  the  trailing  edge  at  t*  = 


Co 


M0-l* 


* 


* 


Figure  7-“  Pressure  distributions  on  wing  undergoing  sudden  angle-of— 
attack  change  at  tf=0. 

For  purposes  of  comparison,  the  growth  with  time  of  the  angle-of— 
attack  indicial  load  coefficient  for  subsonic  flight  is  also  shown  in 
the  part  of  figure  7  entitled  "subsonic."  Since  in  this  case  the 
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lifting  surface  analogue  involves  subsonic  leading  and  trailing 
edges,  the  analysis  requires  the  solution  of  a  boundary-value  problem 
°5.  ^  f®cond  kind-  The  method  of  Eward  (reference  15)  was  used  to 
o  tain  the  results  shown.  It  is  to  be  noted  that  the  expression 


Ap  _  4a 
q  M0 

holds  at  t’=0  for  all  values  of  Mach  number. 

2  Mot' 


Co 


Figure  8.  Indicial  lift— curve  slope  for  Mach  numbers  between 

0  and  1.4  shown  to  time  required  to  travel  12  half-chord  lengths. 

Figure  8  shows  the  variation  of  the  indicial  lift  function 
CLa(tf)  defined  by  the  relation 
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as  a  function  of  Mach  number  and  half-chords  s  traveled  by  the 
airfoil.  The  curve  at  Mb=0  was  first  studied  hy  Wagner  (refer¬ 
ence  22)  and  R.  T.  Jones  (reference  6) .  Since  the  starting  value  is 
it,  CT„(t»)  must  initially  he  infinite.  Immediately  afterward, 

however,  it  assumes  the  value  Jt  and  then  rises  to  the  asymptotic 
value  2  it.  At  a  Mach  number  of  0.4  the  starting  value  of  CLa(t») 
is  10  followed  hy  a  decrease  for  the  time  required  to  travel  approxi¬ 
mately  one-half  chord  length  and  finally  a  steady  rise  takes  place 


to  the  asymptotic  value 


At  Mo=0.8  the  behavior  is  similar 


to  Tne  as  yin  p  ou  tic  vaiuc  — .  ■--«  w 

/l-M02 

The  dashed  portions  of  the  curves  were  determined  from  the  known 
variations  of  the  functions  and  were  not  calculated  explicitly.  The 
asymptotic  values  of  Cia  consistent  with  the  Prandtl-Glauert  correc¬ 
tion  become  so  high,  however,  with  increasing  Mach  number  that  the 
assumptions  of  small  perturbation  theory  are  undoubtedly  invalidated 
near  Mo=l  for  sufficiently  large  values  of  s.  The  initial  portions 
of  the  subsonic  curves  shown  in  the  figure  are,  however,  valid  results 
of  the  theory.  The  nature  of  the  indie ial  lift  function  is  somewhat 
different  at  supersonic  Mach  numbers  in  that  the  beginning  portions  of 
the  curves  are  flat.  The  curves  rise  afterwards,  however,  in  a  finite 
time  to  their  steady-state  value.  From  equations  (4la),  (4lb),  and 
(4lc)  the  expressions  for  CLa(t»)  are  easily  calculated  for  Mo>  1 

and  are  as  follows: 


First  time  interval  0  <  tf  < 


(43a) 


Second  time  interval 


_£o_  <  t*  <  -£o_ 
1+M0  Mq-1 


OLafa*) 


-T— 

it  Mq 


2  +  arc  sin 


Co-Mo* * 


^Mn2-! 


arc  cos 


t»+M0co-t’Mo2 


't»2— (Co— t«Mo) 


(43b) 
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Third  time  interval 


Co 

Mo-1 


<t* 


CLa(tl) 


4 

a/Mo2-1 


(43c) 


Some  of  the  above  results,  along  with  further  developments 
involving  the  entrance  at  supersonic  speed  of  an  unrestrained  ai:n- 
foil  into  a  gust,  have  been  given  in  reference  7- 


Swept— Back  Wings  at  M0=l 

We  turn  now  to  the  special  form  of  the  basic  differential  equa¬ 
tion  for  the  case  M^=l,  As  given  in  table  I,  equation  (E),  the 
velocity  potential  satisfies  Laplace *s  equation 


•  ‘Pyy +  cpzz  =  0 


in  two  dimensions.  The  boundary  conditions  need,  therefore,  to  be 
given  along  strips  normal  to  the  free— stream  direction.  Equation  (7) 
expresses  the  solution  of  the  equation  in  terms  of  two-dimensional 
sources,  sinks,  and  doublets  where  now 


In  rA  =  In  >J (y-ii)2+(z-(;)2 


The  proposed  problem  is  the  determination  of  the  angle— of- 
attack  load  distribution  over  a  swept— back  lifting  plate,  the 
leading  edges  will  be  assumed  straight  lines  while  the  trailing  edge 
will,  for  the  time  being,  be  left  arbitrary.  The  nature  of  the  wing 
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is  thus  indicated  somewhat  arbitrarily  in  figure  9(a)« 


(b)  Plan  form  satisfying  Kutta— Joukowski 
condition. 

Figure  9«—  Swept— back  wings  for  analysis  at  Mq=1. 

Denoting  the  semivertex  angle  by  50  so  that  the  equations  of  the 
leading  edges  are 


y  =  ±  x  tan  50  =  ±  mx 


m 


it  follows  from  equation  (j)  that  since 


dcpu  _  ScPz 

Sz  bz 
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» 

the  velocity  potential  is  given  "by  the  relation 


where 


9(x,y,z) 


Z 

2jc 


A(pn(x,n) 

(y-n)2+z2 


dt] 


Acpo(x,y) 


Avo(x,r),0)  dt) 


mx 


(45) 


It  is  then  possible,  after  integrating  equation  (45)  by  parts  and 
imposing  the  condition  that  A%(x,y)  =  0  at  r\  =  ±  mx,  to  calcu¬ 
late  the  derivative  of  cp  with  respect  to  z  and  thus  obtain  for 
Vo*  induced  vertical  velocity  on  the  wing,  the  expression 


w0 


-1  r ,IGX  Av0(x,ti) 


2  jt 


y— t] 


dT] 


— mx 


(46) 


This  integral  equation  is  to  be  solved  for  AvQ,  the  velocity  w0 
being  assumed  constant  on  the  wing  and  equal  to  — V0a  where  a  is 
angle  of  attack.  The  load  distribution  can  then  be  calculated  from 
this  solution  by  means  of  the  relation 


Ap  _  2  dq>0  _  _2_ 

q  V0  8x  V0  <3x 


Av0(x,r|)  dr) 


(vn 


The  remainder  of  the  analysis  can  best  be  divided  into  two  parts: 
The  first  case  treating  values  of  x  between  0  and  the  second 
dealing  with  the  remaining  values  of  x  on  the  wing. 

Case  I:  0  <,  x  5?  c0.—  Since  the  leading  edges  of  the  wing  are 
of  subsonic  type,  singularities  in  pressure  occur  at  these  edges  so 
that  the  required  solution  of  equation  (46)  is  of  the  form 


Av0(x,y) 


A+By 

Vm^-y2 


(48) 
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Substitution  into  equation  (46)  and  use  of  the  fact  that  Acfo(x,y)  is 
an  even  function  of  y  leads  to  the  result  A=0 ,  B=  —  2wq .  Hence 


I 


« 


and 


Acp0(x,y)  =  -  2w0  v4a2x2-y2 


Ap  _  — 4wpm2x  _  4am2  x 

q  VQ  V m2x2-y2  '/m2x2-y2 


(49) 

(50) 


Case  II:  c0  %  x.-  Let  the  equation  of  the  trailing  edge  be 


y  =  a(x)  or  x  =  a*(y)  (51) 


Using  in  equation  (46)  the  fact  that 


Av0(x,y)  =  -  Av0(x,-y) 


the  expression  for  vq  becomes 


—1  r^(x)  2qAv0dr] 

W°  =  2i  Jo  y2-r|2 


_1_  2TiAv0dq 

25t  va(x)  ^ 


(52) 


If,  on  the  surface  of  the  wing,  Acp0(x,r|)  is  known,  then,  in 
the  wake,  the  discontinuity  in  the  velocity  potential  is  A(P0Ca*(Tl)  J 
since  no  contribution  to  the  Jump  in  potential  is  made  past  the 
trailing  edge.  It  follows  that  if  on  the  wing 


then,  in  the  wake 


0 


Av0(x,q)  =  f(x,T)2) 


(53) 
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* 


Av0(x,t])  = 


^Pq(x,ti)  da*  +  A  dcp0(x,T))  "I 

5x  ^  ^  -*  x=a*  ( r) ) 


(53) 


Jo  Ap[a*(r]),r|]  da* 
2  q  dy 


[a*(T|),ri2] 


(54) 


Substituting  equations  (53)  and  (54)  into  equation  (52)  and  intro¬ 
ducing  the  Kutta  condition  that  loading  at  the  trailing  edge  is  zero, 
-we  get  the  modified  integral  equation 


v 


0 


ffa*(yar),gl3  rV2*2  f(x,ffl) 

a-ffi  ffl  2rt  Ja2  cr-cfi 


dffi 


(55) 


where  the  variables 
The  function 


cr,  hi  now  replace  y2,  t}2, 
f(x,a)  =  2w0 


respectively. 


(56) 


satisfies  equation  (55)  and-  it  remains  now  to  determine  a(x)  so 
that  pressure  is  zero  at  the  trailing  edge.  But  from  equation  (53) 


Acp(x,y)  =  2w0 


yi2— ag(x) 
m2x2-y12 


dyi 


and  thus 


Ap(x,y) 


=  4a 


— k*  —  F(k,t)  +  mE(k,t)  +  -  ft— 

x  V  m 


— k,2m2x2 


dx 


m2x2— y2 


(57) 


where  k’  =  —  and  E(k,\|r)  ,F(k,i|r)  are  incomplete  elliptic  integrals 
with  modulus  k  =  y(l— k'2  and  argument  t=arc  sin  4-  J 1  — 


m2x2 


At  the  trailing  edge  y=a(x)  and 


—  =  4a  [  —  k’  K  +  mE 
q  V  dx 


da 


) 
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where  the  elliptic  integrals  are  now  complete.  If  load  coefficient 
is  set  equal  to  zero,  the  differential  equation 


da  _  m2x  E 
dx  a  K 


(58) 


where 


k  = 


a^ 


follows . 


The  integration  of  equation  (58)  leads  to  the  shape  of  the 
trailing  edge  for  which  the  Kutta  condition  is  satisfied.  Figure  9(h) 
shows  the  plan  form  of  the  wing.  It  can  he  shown  that  the  slope  of 
the  extended  trailing  edge  approaches  the  slope  of  the  leading  edge. 


Figure  10.  Load  distrihution  over  swept— hack  plan  form  at  Mq=1. 

^In  figure  10  the  load  distrihution  is  shown  at  three  spanwise 
stations  for  the  case  when  the  wing  is  cut  off  along  a  line  normal 
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to  the  free— stream  direction.  Over  the  center  section  of  the  wing 
the  Ackeret  type  of  distribution  exists.  The  remaining  sections  have 
discontinuities  in  slope  of  the  loading  at  the  point  where  the  chord 
line  is  cut  by  the  Mach  cone  arising  at  the  trailing  edge  of  the  root 
chord.  This  behavior  of  the  loading  has  been  noted  elsewhere  for 
swept— back  wings  at  higher  Mach  numbers.  (See,  e.g. ,  reference  23.) 

Lift  coefficient  Cp  of  the  wing  is  given  by  the  expression 


where  S  and  b  are,  respectively,  area  and  semispan  of  the  wing  and 
the  first  integral  extends  from  the  leading  edge  to  the  trailing  edge. 
This  equation  may  be  rewritten  as 


2 

Cl  =  g Acp(T.E.,y)  dy 


(59) 


where  Acp(T.E.,y)  is  the  jump  in  potential  at  the  trailing  edge  and 
thus  equal  to  the  circulation  function  r  (y) .  The  following  results 
are  obtained:  ^ 


T  (y)  =  2mV0c0a,  0^y<£A 


r(y)  =  2Voab  Etk^Hl-k^Ftk,*)  j  ,  A^y<b 


(60) 


where  k  and  A  is  the  lateral  distance  to  the  inboard 
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tip  of  the  wing.  (See  fig.  9C13)*) 


In  figure  11  the  value  of 

- £ -  is  plotted  as 

V0c0a  tan  50 


a  function  of  -  ■  ■  ■  ■■■■  1 

Cq  tan  Oq 

for  a  wing  with  semispan 

b=1.325  c0  tan  60  and  aspect 

ratio  AR=4.57  tan  B0.  Results 

of  the  integration  of  equation 

(59)  a**e  shown  in  figure  12 


where 


tan  S0 


is  plotted  as  a 


function  of 


AR 

tan  60 


Figure  11.—  Spanwise  distribution 
of  circulation  for  swept-back 
wing  at  M0=l,  b=1.325  cQ  tan 
and  AR=4.57  tan  50. 


AR/TAN  80 


Figure  12.—  Lift— curve  slope  as  a 
function  of  aspect  ratio  for 
swept— back  wing  at  Mq=1. 


The  methods  presented  here 
can  he  applied  to  the  case  of 
the  swept— back  wing  with  tips 
cut  off  parallel  to  the  free 
stream.  In  this  case  a  Mach 
cone  originates  not  only  at  the 
trailing  edge  of  the  root  chord 
but  also  at  the  intersection  of 
the  leading  and  the  lateral 
edges.  On  the  portion  of  the 
wing  downstream  of  this  Mach 
cone,  the  load  distribution  is 
modified  so  that  an  abrupt  dis¬ 
continuity  exists  at  the  Mach 
cone  and  negative  loading  is 
effective  over  this  part  of 
the  wing.  A  similar  effect  on 
this  type  of  swept-back  wing  has 
also  been  noted  at  higher  Mach 
numbers  in  reference  23* 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif. 
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APPENDIX 

List  of  Important  Symbols 
aQ  speedr- of  sound  in  free  stream 

A  lateral  distance  to  inboard  tip  of  swept  wing  (See  fig.  9) 

AR  aspect  ratio 

b  semi span 

c0  chord  length  (two  dimensions) 

root  chord  (three  dimensions) 

C,  C(0)  load  distribution  factor  introduced  in  equation  (l 6) 

Cl  lift  coefficient 

CLa  lift— curve  slope 

cLa('t0  indicial  lift  function 
Cp  pressure  coefficient 

E  complete  elliptic  integral  of  the  second  kind  with  modulus  k 

Ef ,  Ef0  complete  elliptic  integrals  of  the  second  kind  with  moduli 
a/1—G2 ,  a/  1— 0o2  j  respectively 

Ei,  E2  complete  elliptic  integrals  of  the  second  kind  with  moduli 
ki,  k2,  respectively 


E(k,*) 

P(k,i|r) 


G  parameter  defined  in  equation  (29) 

V  Cq9q  • 

x  c0 
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complete  elliptic  integral  of  the  first  kind  with  modulus  k 


Ki,  K2  complete  elliptic  integrals  of  the  first  kind  with  moduli 
ki>  k2,  respectively 

m  tan  60 

Mq  free— stream  Mach  number 

p  local  static  pressure 

p0  free— stream  static  pressure 

q  free— stream  dynamic  pressure  Qjp0V 02^) 

7(x-|)2+p2(z-^)2 


(x-e)2+32[(y^)2+(z_^)2] 

(x^)2-p2[(y_r1)2+(&_f)2] 


distance  traveled  in  half— chords 


area  of  wing 


u,v,w  perturbation  velocity  components  parallel  to  x,y,z  axes, 
respectively 

Au0,  Aw0  jump  in  value  of  u,  w  at  the  z=0  plane 
VQ  velocity  of  free  stream 

x,y,  z  Cartesian  coordinates 

xf,zf,tf  coordinates  introduced  in  equation  (2) 


angle  of  attack  in  radians^ 

1 1-Mo2  | 

load  coefficient  ^C  PZ^Pu) 
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A  semivertex  angle  of  yawed  triangular  wing 

8  angle  "between  lifting  element  and  x  axis 

&o>  &i  angles  between  leading  edges  of  yawed  triangle  and  x  axis 


Q>Qo,Q-i.  P  tan  6,  p  tan  S0,  p  tan 
A  sideslip  angle  of  yawed  triangle 


4 

Po 

T 

<P 

05 


Mach  angle  (^arc  sin 

Cartesian  coordinates 

free— stream  density 

region  of  integration  in  equation  (10) 

perturbation  velocity,  potential 

pZ 

sign  denoting  "finite  part"  of  integral 

Subscripts 


subscript  denoting  value  of  variable  on  upper  surface 
of  ving 

subscript  denoting  value  of  variable  on  lower  surface 
of  wing 
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